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1 Introduction 



According to the pioneering paper |J by Haag on the BCS theory of super- 
conductivity p, the general meaning of Cooper pair condensates || can be 
understood as a "variable at infinity" 0] which is based on the cluster prop- 
erty valid in any thermodynamic pure phases described by factor states^. 
Along this line, Sewell || takes, as a general characterization of the super- 
conducting BCS states (-)bcs with condensed Cooper pairs, the off-diagonal 
long range order (ODLRO): 

\mxi + §m*i - f )^(* 2 - + §)w 

-vix^^ix^)] -+ 0. (1) 

IX1-X2I-00 



1 A factor state means a state such that the (GNS) representation of the algebra of 
physical quantities corresponding to it has a trivial centre. 
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Here ip(x) denotes the second-quantized non-relativistic electron field obeying 
the canonical anticommutation relation 



{1>(x), ^(y)} = 6 3 (x-y), (2) 

and the so-called "macroscopic wave function" £) = (ip(X + §)^(A" — 

I)) bcs of a Cooper pair should be non- vanishing in the limit of spatial infinity 
\X\ — > oo. According to this formulation, he succeeded in giving a general 
proof of the validity of M eissner effect (B = inside of superconductor) . 

Inspired by these attempts to understand the essential features of super- 
conductivity in a model-independent way, I try here to present a simple and 
general derivation of Josephson formulae?] || which describe the tunnelling 
currents (dc and/or ac) caused by the phase differences between two super- 
conductors separated by a thin barrier of the insulator. In contrast to the 
traditional derivations based on the tunneling Hamiltonian ||, we see that 
they are just a simple and direct consequence of the above ODLRO and of 
the fact that the energy contained in one side of the junction gives a non- 
trivial response to the global gauge transformation caused by the presence 
of phase difference. 

2 Simple Derivation of Josephson Formulae 

Although a completely model-independent approach is desirable, the present 
discussion requires the postulate of the standard BCS Hamiltonian [fj, [l[ 
arising effectively from the electron-phonon coupling: 

H BCS (A;w) = I dx[^-V^(x) ■ Vip(x) - ^(x)ij(x)} 
J a 2m 

+ — dx dy dz du w(x,y, z 1 u)^(x)tjj ji (y)ijj(u)ip(z), (3) 

\M J A J A J A J A 

where m is the mass of an electron, fi the chemical potential of electrons 
in the bulk superconductor occupying a macroscopically extended spatial 
region A with volume |A|. To exhibit the essence, we use here a simple 
idealized picture of weakly coupled superconductors placed in two spatial 
regions A 1; A 2 (c A) separated by a Josephson junction regarded ideally 
as a phase boundary W = dAi = dA 2 . In view of the wide applicability 
of BCS model (at least in non high T c cases), the possible differences in the 

2 While this issue is treated in || , the arguments there are not clear or explicit enough. 
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properties of superconductors in Ai and in A 2 are supposed to be absorbed in 
the different choices of the potential functions Wai (x, y, z, u) = w(x,y, z,u) 
for x, y,z,u G Ai and wa 2 (x, y, z, u) = w(x, y, z, u) for x, y, z,u G A 2 . Aside 
from this freedom, we understand that the dynamics of superconductors are 
described universally by the above H BC s an d that the differences of the 
realized thermodynamic phases are all reduced to those in the choice of states 
(•)• 

Here a remark need be added on a subtle point: To give a precise meaning 
to the term "thermodynamic phases" , one should consider the theory in the 
thermodynamic limit with volume tending to infinity. Namely, the sizes of 
regions Ai, A 2 (of course, finite on the macroscopic scale) should be treated 
as "infinitely large" according to the scale of microscopic interactions, and 
the location of the junction should be supposed to be at "spatial infinity" 
far away from (the centre of) Ai in this scale. Although such expressions as 
this may sound quite naive and vague lacking in the mathematical rigour, it 
is possible to give a mathematically precise meaning to it in a non-standard 
analytic framework as will be discussed briefly in the next section. 

The next essential ingredient is the very definition of the tunneling electric 
current flowing through the barrier. While our system does not have a locally 
conserved electric current owing to the presence of non-local coupling in 
Eq.(^]), the conservation of electric current is still meaningful in the following 
sense: We define the electric charge Q(A) in a spatial region A by 

Q(A) = -|e| / dx ip\x)ip{x), (4) 

J A 

which is conserved in the sense of 

[H BCS (A;w), Q(A)] = 0, (5) 

and which generates the electric global {/(l)-gauge transformation: 

[iQ(A), V(x)] = i\e\il>{x), [iQ(A), ^(x)} = -i\e\^(x), (6) 

where e = — \e\ is the unit of electric charge. In the situation with A = 
Ai U A 2 (UW 7 ), the tunneling current J between the two regions Ai and A 2 
is defined by 

J = ^(QiA^Bcs = ([iH BCS (A), Q(A 1 )]) BCS . (7) 

Note that J cannot be non-vanishing without the presence of the region A 2 
outside of Ai in view of the simple equality Eq.(||) applied to A = Ai. (In 
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du 



a sense Eq. (|7|) can be viewed as the basis of heuristic expression dN / dt = 
dH/d6 discussed in the number-phase picture of Ginzburg-Landau theory.) 

In view of the local (anti-)commutativity following from Eq.(0), the right- 
hand side of (0) can be reduced to 

J = ([iH 12 , Q(Ax)])bcs = -<[i'Q(A a ), H 12 ]) BCS . (8) 

Here, H\ 2 appears in the following decomposition of Hbcs(A; w) correspond- 
ing to that of the spatial region A into Ai and A2 (separated by W): 

H BCS {A; w) = H BCS {Ai] + H bcs{M] jj^w) + H l2 , (9) 

according to which it can be written explicitly as 

H± 2 — ttt ^ f dx f dy f dz f du w(x,y, z,u)^ {x)%p ] (y)ip(u)ip(z). 

1 \{i,j,k,l}={l,2} Ai Aj Afe Ai 

(10) 

Then the commutator in Eq.(H) is calculated as 

- ([iQ(Ai), H 12 \) BC s = ^Mttt ^2 J dx I dy J dz 

1 Iy \{i,j,k,i}={m Ai Aj Ak Al 

x w(x,y,z,u)(-Sn - 5 jX + hi + 5n)(^ ] {x)ip ] {y)i){u)4){z)) BCS - (11) 

To give the precise meaning of the "macroscopic wave function" \I/(X, £) in 
ODLRO (JJ), we should now consider the thermodynamic limit of A, Ai, A2 — ► 
00, in which we assume the ratio |Ai|/|A2| is kept fixed. When we evaluate 
the right-hand side of (^) in this situation, the contributions of terms coming 
from the region A 2 can be replaced by their expectation value in the BCS 
state (-)bcs because of the cluster property justified by the understanding 
that the boundary W is infinitely far away in this limit, A 1; A 2 — > 00 with 
|Ai|/|A 2 | kept fixed: 

■ --**{x®)4#{xf) ■ ■■^*(x%)) BCS 

- &*(x?) ■ ■■^*(x^)) BC s^ # (x { 1 2) ) ■ ■■^*(x^))bcs\ 

— > 0, (12) 

Ai,A 2 — >oo, \Ai\/\A 2 \:fixed 

where x^\--- ,x$ G Ai, xf\--- G A 2 and ip# = ip or ift. Since 

the expectation values of odd powers of fermionic operators ip, if)' van- 
ishesQ, the only contributions surviving this limit to the above integrand 

3 This holds under the moderate assumption of spatial homogeneity at infinity [Q 
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XA i (x)xA j (y)XA k (z)xA l (u)(-8u - 6ji + 5 k i + 5 n )(tfj ] \x)^ '{y)tp{u)tfj(z)) BCS 
come from the cases with i — j — l,k — I — 2ovi— j — 2, k — I — 1: 



x(x,y, z,u)(-5n - Sjt + S k i + 5n)(ip\x)^(y)ip(u)ilj(z))Bcs 

{i,j,k,l}={!,2} 

— ► X(x,y,z,u)(-6 il - 6^ + 6111 + S a ) 

Ai,A2 — >oo, Ai \/\A2\- fixed z — * 

i=jj^k=l 

x {^{x)^{y)) BC s^{u)^{z)) B cs , (13) 

where x{^,y,z,u) = XaX x ) XA j (y)XA k {z)XA l {u) with xa being the indicator 
function of a set A. If we assume the almost spatial homogeneity^ in each of 
regions Ai, A 2 and if the potential w(x, y, z, u) between electron pairs can be 
assumed to be a real quantity, we can put 

i , f \ ,f w f l^ik^ 1 (x,y G Ai), M 

M x Mv))bos*[ l^^je^j; (14) 

from which we obtain the desired result for the dc- Josephson current by 
picking out the phase factors from the integrands: 

J ~ constant x ( e i(9l ^ 2) - e ~ i{dl ' 62) ) oc sin(A0), (15) 

where A6 = 61 — 62 is the phase difference of Cooper pairs across the junction. 

Once this is obtained, it is straightforward to derive the formula for ac- 
Josephson current applied to the situation with voltage gap V across the 
junction by replacing the above A6 with A6 + 2eVt. While the postulated 
BCS Hamiltonian (§) with non-local interactions is not fully compatible with 
a local gauge invariant coupling of the system with spatially varying external 
electromagnetic field, it still allows a coupling with V in a way invariant under 
local gauge transformations with only temporal dependence. As argued in 
0, this gauge freedom allows us to treat the situation with voltage gap V 
simply by performing a time-dependent local gauge transformation: 

A» = (</>, A=0)^(<f> + ^(Vt) = V, 1=0); (16) 
ip(x,t) — ► e ieVt i/;(x,t), (17) 

which causes the above shift A6 — ► A6 + 2eVt of the phase difference. Thus 
we have the formula J ac oc sin(A# + 2eVt) for ac- Josephson current. 

4 This can be formulated without difficulty in the non-standard analytic framework 
mentioned in Sect. 3 in such a way that suitable combinations of quantities involving the 
deviations from the homogeneity arc infinitesimal. 



5 



In deriving Eq.(|T5D we need the spatial homogeneity to extract the factor 
of phase difference. Although the above kind of qualitative discussion does 
not allow us to determine the precise coefficient, we can extract the contri- 
bution to the energy density of the phase difference at the boundary located 
infinitely far away in the similar way to the above, which gives us 



( Hl flf cs — > constant x cos(A#). (18) 

|A| Ai,A 2 — >oo, |Ai|/|A 2 ]:/iised 

If the coefficient of cos(A#) is assured to be of negative sign, this guar- 
antees the self-consistency of the postulate that the phase of Cooper pair 
condensates in a superconductor should be spatially homogeneous in favour 
of A6 = in the absence of such a constraint to maintain the phase difference 
as the barrier. To verify this consistency problem in a more satisfactory way, 
it seems necessary to confront a challenging problem of how one can justify 
the notion of a point-like order parameter ty(x) of Cooper pairs, which ap- 
pears in the Ginzburg-Landau phenomenological approach ||, and which is 
crucial for discussing the Type II superconductivity involving in an essential 
way the spatial inhomogeneity and the local gauge invariance problem.pl 



3 Discussion 

Although we refrain from a systematic explanation, we comment here how the 
method of non-standard analysis || can be useful in describing the situation 
with infinitely large regions A 1; A 2 separated by a boundary W at infinity. 
What is important is that it allows us to treat both the finite volume theory 
and the infinite volume one simultaneously without disconnecting the two 
approaches. Before introducing the distinctions among finite, infinity and 
infinitesimal (at the syntactic level of internal objects), everything looks as if 
we were in the finite volume theory, but, once such distinctions are introduced 
(by the interpretation in a non-standard model) by regarding |Ai| and |A 2 | 
as infinite numbers (whose ratio is kept finite), the infinite volume theory is 
seen to be contained in the former, being extracted through the procedure 
of taking finite parts of quantities which throws away all the infinitesimals 
such as l/|Ai|. In the situations with only one thermodynamic phase, this 
kind of treament does not make much difference from the usual one, because 
the procedure of extracting finite parts is more or less equivalent to taking 
the thermodynamic limit. In the present situation with two infinitely large 

5 This sort of situations may be the very places where the conceptual advantages of 
non-standard analytic formulation should be exhibited. 



6 



regions Ai and A2, however, we have still "another world" in A2 beyond 
the infinitely distant boundary W of infinitely extended A^. Whereas this 
situation seems difficult to be accommodated in the usual formulation, it can 
be described without difficulty in the framework mentioned here, where all 
the infinities and infinitesimals are fully legitimate quantities. Moreover, all 
such limiting or approximate relations as Eqs. ([!]), fljjD, (pT3[), ((Hp, (|i"5|), (p~8|) 
are replaced by simple algebraic equivalence relations modulo infinitesimals, 
in which one of its conceptual advantages can be found. 
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